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Lithium-rich field giants in the Sloan Digital Sky Survey 
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ABSTRACT 

We present a search for post-main-sequence field stars in the Galaxy with atypically large 
lithium abundances. Using moderate-resolution spectra taken as part of the Sloan Digital 
Sky Survey, along with high-resolution followup spectroscopy from the Hobby-Eberly Tele- 
scope, we identify 23 post-turnoff stars with loge(Li) greater than 1.95, including 14 with 
loge(Li) > 2.3 and 8 with loge(Li) > 3.0, well above the low level expected for evolved 
stars. Comparison with theoretical isochrones indicates that some of our Li-rich stars are af- 
filiated with the upper red giant branch, the asymptotic giant branch and the red clump rather 
than the RGB bump, which is a challenge to existing models of Li production in evolved stars. 
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1 INTRODUCTION 

The abundance of lithium is extremely sensitive to the ambient con- 
ditions inside stars, since its cross section to proton capture is quite 
large at temperatures above 2.6 x 10 6 K. Li is readily destroyed in 
the interiors of low-mass main-sequence stars, and the Li in the 
convective envelope is progressively depleted as envelope mate- 
rial is cycled past the innermost, hottest region of the convective 
zone. Surface Li abundance is further deplet ed after star s leave the 
main sequence and undergo first dredge-up (Iben 1968), in which 
the convective envelope deepens, mixing fusion-processed material 
from the stellar interior into the envelope and exposing envelope 
material to even higher temperatures. 

Because of t his, it was quite surprising when 

IWallerstein & Snedenl J 19820 discovered a Li-rich K giant 
serendipitously while investigating barium-ri ch field stars . More 
systematic searche s for Li-rich stars (e.g., iBrown et all Il989l: 
Pilach owskil 1 19861 : IPilachowski et alj 1 19881 : IPilachowski et ail 
that they are rare in the Solar neighborhood 
in globular clusters. The study of 
20000 found that nearby Solar- 
1M ) 



20001) found that they 



and extremely unusual 
ICharbonnel & Balachandranl 

metallicity Li-rich stars are either low-mass stars (M 
near the "bump" in the red giant branch (RGB) luminosity function 
or intermediate-mass stars (M ~ 3Mq) low on the asymptotic 
giant branch (AGB). In both cases, the hydrogen-burning shell 
in the stars in question has recently encountered, and destroyed, 
the discontinuity in mean molecular weight (the "/i-barrier") 
left behind at the end of first dredge-up. In low-mass red giants, 
the destruction of the /i-barrier causes a brief loop in stellar 
evolution, observed as the RGB bump, and permits deep mixing, 
a nonconvective process that steadily cycles material between the 
photosphere and the radiative zone, causing surface abundances of 
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carb on and nitrogen to e volve constantly as stars ascend the RGB 
(see iMartell et al.1l2008l for a more thorough explanation of deep 
mixing). In higher-mass stars, the destruction of the /i-barrier oc- 
curs at second dredge-up, low on the asymptotic giant branch. The 
fact th at the Li-rich stars identified bv lCharbonnel & Balachandranl 
(2000) were located near the destruction of the /i-barrier led those 
authors to speculate that the onset of deep mixing produces a brief 
episode of Li production via the Cameron-Fowler mechanism 
dCameron & Fowlerlll97ll) . In this process, 7 Be produced by the 
proton-proton chain is transported to a cooler zone in the star, 
where it captures an electron and decays into 7 Li, which is stable 
against proton capture if its surroundings are cool enough. A phase 
in which the deep-mixing mass transport speed was briefly high 
enough to allow this process would cause surface Li abundances to 
rise sharply and then decay. 

However, the existence of Li-rich stars in other regions of 
the He rtzsprung-Russell diagram challenges this model: iKraft et al.1 
(1999) found one bright red giant well above the RGB bump in the 
globular cluster M3, where deep mixing is well known to operate 
(Smith 2002: lAngelou et alj|201ll) . More recently, six independent 
studies in the past three years have reported a total of 47 new Li- 
rich post-main-sequence stars. These new Li-rich stars fill a broad 
range of parameter space: low-mass (M < 2M^) RGB sta rs in 
the thick disk with -1 < [Fc/H] < dMonaco et all 1201 ll) : an 
inter mediate-mass (M ~ 4Mq), Solar-metallicity star in the thin 
disk ^Monaco et ail 1201 ll) : l ow-mass RGB stars with [Fe/H] ~ 
—0.9 in the Galactic bulge dGonzalez et"aH [2009: Lebzelt er et al.1 
120121) : and low-mass RGB stars in Local Group dw arf galaxies, 
two of which have [Fe/H] ~ -2.7 dKirbv et alj2012t) . 

T wo of the recent studies ( Kuma r~et alj|20 1 ll and lRuchti et al.l 
l201lh take the approach of selecting a well-defined subset from 
an existing catalog (Hipparcos and RAVE, respectively) for high- 
resolution spectroscopic followup and abundance analysis, a pro- 
cedure that we also adopt. In this paper, we present a search for Li- 



2 S. L. Martell & M. D. Shetrone 





Figure 1. Distribution of T e g and log(g) for the 8535 stars in our initial 
data set. 



rich field stars in Sloan Digital Sky Survey (SDSS) spectroscopy. 
By surveying stars across a wide range of metallicity and surface 
gravity, we are able to take a broad view of the question of Li-rich 
post-main-sequence stars. This paper describes the 23 previously 
unknown Li-rich stars we identified from this original data set, and 
explores whether the various proposed mechanisms for Li produc- 
tion in evolved stars are sufficient to explain this group of Li-rich 
stars. 



2 THE OBSERVATIONAL DATA SET 

In order to explore the physical conditions that produce large 
amounts of photospheric lithium in some evolved stars, we selected 
a sample of post-main-sequence stars from the SDSS Data Release 
7. The 8535 stars in our initial data set were required to have high- 
quality spectra and relatively low errors on stellar param eters deter- 
mined by the automated SSPP pipeline iLee et all2o"o83) . and were 
allowed to cover a broad range of metallicity, effective temperature 
and surface gravity. 



2.1 SDSS spectroscopy 

The imaging component of the origina l Sloan Digital Sky Sur- 
vey and its extensions jFukugita et al. 1 1 19961: Gunn ct apl 19981; 
York et al.l2000t|Pier et al.l2003LlGunn et al 2006; Sto ughton et ail 



2002 | ;lAbazaiian et alBoolll2004l2005l;lAdelman-McCarthv et al.1 
20061 12007L 120081 ; lAbazaiian et alj|2009l ; lAihara et alj|201lh "h> 
eludes uqriz phot ometry for several hundred million stars. SEGUE 
dYannv et alj|2009h , one of three sub-surveys that together formed 
SDSS-II, extended that imaging footprint by approximately 3500 
deg 2 , and also obtained R ~ 2000 spectroscopy for approximately 
240, 000 stars over a wavelength range of 3800 - 9200A. Ad- 
ditional spectroscopy with similar ch aracteristics has bee n taken 
as part of the SDSS Legacy project ( Abazai ian et al. I I2009I) and in 
preparation for the SDSS-III Multi-object APO Ra dial Velocity Ex - 
oplanet Large-area Survey (MARVELS) project dGe et alj|2008h . 
To facilitate the development and calibration of the SEGUE Stel- 
lar Parameter Pipeline (SS PP; iLee et al]|2008at iLee et al.ll2008bl ; 
lAllende Prieto et al.1 I2008T) . SEGUE included spectra for a col- 
le ction of Galactic globular and open clusters. As described 
in ILee et aD 42008a), SSPP produces estimates of T c h, log(g), 



Figure 2. Metallicity histogram for our initial data set (heavy solid line), 
with histograms for the three data sources overplotted with a dotted line 
(SEGUE), light solid line (Legacy) and dashed line (MARVELS). 



[Fe/H], and radial velocity, along with the equivalent widths and/or 
line indices for 85 atomic and molecular absorption lines, by pro- 
cessing the calibrated spectra g enerated by the standard SDSS spec- 
troscopic reduction pipeline dStoughton etail |2002|) . Additional 
validations and details of recent updates to SSPP are discussed in 
ISmolinskietal.U201lh . 



2.2 Target selection 

We selected our sample from the seventh SDSS data release (DR7), 
using the online Catalog Archive Server (CAs£] to identify stars 
with a mean signal-to-noise ratio per pixel (SNR) of at least 40, 
T c ft < 6000K, reasonable errors on [Fe/H] metallicity (o"[f /h] < 
0.3) and surface gravity (o"i og ( g ) < 0.5), and surface gravity no 
higher than 1 .4 dex greater than the surface gravity at the "bump" 
in the RGB luminosity function. Previous studies have found that 
stars with unexpectedly high lithium abundances tend to be near 
the RGB bump, so we designed this gravity selection to include 
stars at evolutionary phases from shortly before the RGB bump 
to the tip of the RGB. Empiricall y, the luminosity of th e RGB 
bump is a function of metallicity: iFusi Pecci et al. (1990) report 
[Fe/H] versus A4"^ ump for a selection of Galactic globular clusters, 
and we converted those My values to lo g(g) using metallicity - 
matched 12 Gyr Yale-Yonsei isochrones (Demarque et al.ll2004i) . 
then fit a least-squares line to derive log(gJ 



bump 



= 2.23235 + 



(0.188018 x [Fe/H]). This initial data set contained 8535 stars: 
7667 from SEGUE, 471 from Legacy and 397 from MARVELS. 
Figure 1 shows T c ff versus log(g) for all 8535 stars, and Figure 2 
is a histogram of the metallicities. In Figure 2, the heavy solid his- 
togram shows the metallicity distribution for the full data set, with 
the histogram for just the SEGUE stars shown as a dotted line, the 
histogram for the Legacy stars drawn with a light solid line, and the 
histogram for the MARVELS stars shown with a dashed line. 



3 LOW-RESOLUTION ANALYSIS 

In order to identify stars with potentially high lithium abundance 
based on low-resolution spectra, we designed a spectral index, 
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Figure 3. Lithium index S(Li) versus T e fj for the initial data set (upper 
panel); S(Li) versus [Fc/H] (lower panel). In both panels, stars selected 
for the secondary data set are shown as open triangles. 




Figure 4. Typical spectra from our initial SDSS low-resolution data set. 
The lower two spectra appear to have strong absorption in the 6708A Li 
resonance line, while the top spectrum does not. 



S(Li), that measures the magnitude difference between the inte- 
grated flux in the 6708A resonance line and the integrated flux in a 
small region of nearby spectrum. It is defined as 

I\dX 

S(Li) = -2.5 xlog j6 ™> 

J6730 l * dX 

Typically spectral indices are defined for use in small regions of 
parameter space, and use one or more nearby region(s) of spectrum 
independent of the property that drives variations in the feature of 
interest as a comparison band; however, our initial data set con- 
tains stars across such a wide range of T e g, log(g) and [Fe/H] that 
there are no nearby regions of spectrum that are the same for all of 
our stars. This fact introduces some scatter to our S(Li) measure- 
ments, requiring extra care in selecting likely Li-enhanced stars for 
followup observations. 

The upper panel of Figure 3 shows S(Li) versus T e g for the 
initial data set, and the lower panel shows S(Li) versus [Fe/H] 
for the same stars. In each plane, S(Li) has a broad distribution 
at a fixed value of the independent parameter (T e ff or [Fe/H]) and 
a slight dependence on that parameter. To identify outliers in the 
S(Li)-T c g and S(Li)-[Fe/H] distributions, we select all stars with 
S(Li) > 0.295 - jSfJr and S(Li) > 0.18 + as our sec- 

ondary data set. There are 162 stars meeting both criteria, and these 
are plotted as open triangles in Figure 3, while all other stars from 
the initial data set are represented by small filled circles. 

From by-eye examination each of the spectra in this secondary 
data set of 162 stars, we found that our selection criteria were gen- 
erous: fewer than half of the stars selected based on S(Li) - T e ff 
- [Fe/H] position had 6708A Li lines that appeared legitimately 
strong in the SDSS spectra. Figure 4 shows spectra for two of 
the stars we selected for high-resolution followup (the lower two 
spectra), as well as one star with little apparent absorption in the 
6708A Li line (top spectrum). We generated low-resolution syn- 
thetic spectra based on SSPP stellar parameters to roughly estimate 
lithium abundance for those potentially lithium-enhanced stars, and 
selected 36 of them for high-resolution followup observations, with 
the goal of determining more precise stellar parameters and e(Li) 
abundances. Table 1 lists SDSS ID number, data source (SEGUE, 
Legacy or MARVELS), plug-plate number, modified Julian date of 
observation, fiber ID number, right ascension, declination, g mag- 



nitude (V magnitude for the MARVELS stars) and S(Li) measure- 
ment for the 36 stars initially flagged for followup spectroscopy. 
MARVELS stars are too bright to be included in SEGUE photom- 
etry, so their apparent V magnitudes were taken from the Vizier 
catalog-search service Q. The plate, MJD and fiberlD information 
is given because it is a simple way for finding individual stars in 
the CAS. 



4 HIGH-RESOLUTION ANALYSIS 

Followup spectra were obtained for 34 of these 36 stars using the 
High-Resolution Spectrograph (HRS) o n the Hob by-Eberly Tele- 
scope (HET) at McDonald Observatory jTullll 1998h . All but one of 
these 34 stars we re first observed as fi ller time (P4) during normal 
queue observing dShetrone et ai]2 007 ) during periods 2010T2 and 
201 IT 1. This filler time allowed us to take an initial short snapshot 
spectrum to confirm the stength of the Li line; 7 were found to not 
be Li rich and 27 were confirmed Li rich. For 23 of these Li rich 
stars we were able to follow-up with additional spectra to get the 
SNR sufficient for determination of stellar parameters and abun- 
dances. One of the 36 stars (SDSS J083 1+5402) had already been 
abserved as part of the Chemical Abundances of Stars in the Halo 
(CASH) project dFrebel et al.l2008l) . 



4.1 Spectroscopic Observations and Reduction 

HRS was configured to HRS_15k.central_316g5936.2as_0skyJS0.GC0_2x5 
during our HET observations to achieve R= 18,000 spectra covering 
4233A to 5900A on the blue chip and 6000A to 7850A on the red 
chip. Exposure times per visit ranged from 300 seconds to 1800 
seconds based on the observing conditions. Because of the nature 
of the HET and the generally poor conditions of filler time, some 
targets had multiple visits taken over a number of days to build up 
enough S/N to cleanly detect the Li feature. 

The spectra were reduced with IRAH3 ECHELLE scripts. 



2 http://vizier.cfa.harvard.edu/viz-bin/VizieR 

3 IRAF (Image Reduction and Analysis Facility) is distributed by the Na- 
tional Optical Astronomy Observatories, which are operated by the Associ- 
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Table 1. Star identifiers, position, photometry and lithium index for the 36 stars selected for followup spectroscopy. Apparent V magnitudes are listed for 
MARVELS stars. 



SDSS ID 


Source 


Plate 


MJD 


FiberlD 


a 


S 


90 


S(Li) 








Confirmed Li-rich with HET 








SDSS J2206+4531 


SEGUE 


2556 


54000 


640 


22:06:19.75 


45:31:57.19 


15.732 


0.209 


SDSS J2353+5728 


SEGUE 


2377 


53991 


337 


23:53:29.76 


57:28:18.41 


12.960 


0.209 


SDSS J0808-0815 


SEGUE 


2806 


54425 


256 


08:08:34.33 


-08:15:47.30 


14.598 


0.198 


SDSS J0301+7159 


MARVELS 


2843 


54453 


634 


03:01:33.77 


71:59:11.16 


10.570 


0.197 


SDSS J2019+6012 


SEGUE 


2554 


54263 


600 


20:19:32.44 


60:12:46.32 


15.270 


0.212 


SDSS J0652+4052 


MARVELS 


2847 


54452 


557 


06:52:37.80 


40:52:01.44 


11.170 


0.212 


SDSS J0245+7102 


MARVELS 


2843 


54453 


168 


02:45:20.83 


71:02:02.63 


11.760 


0.184 


SDSS J0632+2604 


SEGUE 


2678 


54173 


188 


06:32:30.90 


26:04:37.46 


15.152 


0.234 


SDSS J2356+5633 


SEGUE 


2377 


53991 


200 


23:56:52.22 


56:33:38.83 


13.080 


0.210 


SDSS J0535+0514 


MARVELS 


2841 


54451 


299 


05:35:43.75 


05:14:22.62 


12.310 


0.191 


SDSS J2200+4559 


SEGUE 


2556 


54000 


540 


22:00:51.22 


45:59:40.33 


15.453 


0.229 


SDSS J0304+3823 


SEGUE 


2441 


54065 


453 


03:04:37.40 


38:23:46.12 


15.357 


0.181 


SDSS J1901+3808 


SEGUE 


2536 


53883 


306 


19:01:49.00 


38:08:51.00 


14.29 


0.214 


SDSS J0720+3036 


SEGUE 


2677 


54180 


255 


07:20:42.79 


30:36:51.49 


16.667 


0.209 


SDSS J1909+3837 


SEGUE 


2536 


53883 


110 


19:09:29.25 


38:37:46.23 


15.764 


0.214 


SDSS Jl 105+2850 


SEGUE 


2870 


54534 


599 


11:05:36.91 


28:50:08.09 


17.006 


0.185 


SDSS J0654+4200 


MARVELS 


2847 


54552 


614 


06:54:58.00 


42:00:27.00 


12.51 


0.220 


SDSS J1607+0447 


SEGUE 


2178 


54629 


378 


16:07:09.23 


04:47:12.73 


16.139 


0.173 


SDSS J1310-0012 


SEGUE 


2901 


54652 


191 


13:10:37.22 


-00:12:44.40 


15.252 


0.185 


SDSS J0936+2935 


SEGUE 


2889 


54530 


164 


09:36:27.44 


29:35:35.80 


16.240 


0.185 


SDSS J1522+0655 


SEGUE 


2902 


54629 


33 


15:22:11.86 


06:55:55.22 


15.912 


0.165 


SDSS J083 1+5402 


SEGUE 


2316 


53757 


395 


08:31:55.28 


54:02:45.55 


15.430 


0.186 


SDSS J1432+0814 


SEGUE 


2933 


54617 


107 


14:32:07.15 


08:14:06.15 


17.047 


0.221 



Confirmed Li-rich with HET but not enough SNR for full analysis 



SDSS J0515+1558 


SEGUE 


2668 


54084 


214 


05:15:23.17 


15:58:55.38 


14.555 


0.216 


SDSS J2048+5603 


SEGUE 


2555 


54265 


58 


20:48:46.82 


56:03:53.65 


13.423 


0.200 


SDSS J21 19-0734 


SEGUE 


2305 


54414 


230 


21:19:47.65 


-07:34:26.89 


14.526 


0.167 


SDSS J1214+5027 


SEGUE 


2919 


54537 


630 


12:14:54.88 


50:27:00.99 


18.205 


0.177 



Confirmed not Li-rich with HET 



SDSS J1421+0136 


Legacy 


534 


51997 


162 


14:21:28.26 


01:36:28.35 


17.077 


0.200 


SDSS J1213-0147 


Legacy 


333 


52313 


331 


12:13:53.32 


-01:47:14.93 


16.558 


0.180 


SDSS J1227-0026 


SEGUE 


2558 


54140 


115 


12:27:12.72 


-00:26:58.81 


16.990 


0.177 


SDSS J0144-1005 


SEGUE 


2850 


54461 


12 


01:44:49.17 


-10:05:31.09 


16.738 


0.171 


SDSS J0352-0540 


SEGUE 


2051 


53738 


357 


03:52:43.59 


-05:40:03.32 


16.487 


0.168 


SDSS J1725+0729 


SEGUE 


2797 


54616 


304 


17:25:04.06 


07:29:45.03 


14.807 


0.168 


SDSS J1849+1931 


SEGUE 


2812 


54633 


393 


18:49:52.58 


19:31:04.20 


14.980 


0.168 



Not observed with HET 

SDSS J1824+7911 SEGUE 2802 54326 372 18:24:22.51 79:11:24.47 16.548 0.201 

SDSS J2021+1257 SEGUE 2248 53558 291 20:21:25.73 12:57:51.33 15.263 0.198 



The standard IRAF scripts for overscan removal, bias subtrac- 
tion, flat fielding and scattered light removal were employed. 
For the HRS flat field we masked out the Li I and Na D re- 
gions because the HET HRS flat field lamp suffered from inter- 
mittent faint emission lines which would otherwise compromise 
our analysis. Equivalent widths (EW) were measured using the 
IRAF task splot and forcing the Gaussian fit to have fixed FWHM 
for those lines below 150 mA while allowing the strong lines 
to have a variable FWHM. We also fit multiple Gaussian pro- 
files to any nearby lines to take blends into account. Metallic- 
ity and a-element abundances were determined through equiva- 
lent width analysis; oscillator strengths ad o pted for the EW analy- 
sis were taken from lFulbright et al.l(l2006l) . iFulbright et ail ( 120071) , 



ation of Universities for Research in Astronomy, Inc., under contract with 
the National Science Foundation. 



Ijohnson & Pilachowskil bold) , and iKoch & McWilliaml (2008). 
Lithium abundance was determined by spectral synthesis of the 
6708 A and 6 104 A resonance lines. We also determined carbon 
abundance and 12 C/ 13 C ratios, where possible, by synthesis of 
the 4297A CH feature. The CH linelist for the carbon analysis 
was taken from Plez (priv. comm.); all other lines were taken from 
a version of the iKurucz & BelH £l995) linelist modified to fit the 
strongest lines to the Solar spectr um. Model atmosphere s were 
taken from the Osmarcs 2005 grid (Gusta fsson et al.l l2008). inter- 
polating between the grid points. These spherical models cover 
metallicity ranges from +1.0 to -3.0 with [a/Fe] = +0.4 for 
[Fe/H] <= -1.5, [a/Fe] = +0.3 for [Fe/H] = -0.75, 
[a/Fe] = +0.2 for [Fe/H] = -0.5, [a/Fe] = +0.1 for 
[Fe/H] = -0.25 and [a/Fe] = +0. for [Fe/H] >= 0.0). 

The abundance calculations were performed using a modified 
version of the 2010 versi on of the LTE line analysis and spectrum 
synthesis code MOOG dSneden|[l973l) . This code has some sig- 
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Figure 5. The 23 stars for which we derive high-resolution stellar parame- 
ters and abundances, in the a - Fe plane. Dashed lines mark the borders of 
our rough age estimations. 



nificant differences from the earlier versions including using the 
predicted abundance vs. reduced EW for determining the microtur- 
bulence. Initial estimates of the stellar parameters were taken from 
the SSPP pipeline. From this starting point the temperature was 
adjusted until the there was no significant slope in the Fe I abun- 
dances with respect to line excitation potential. The stellar gravity 
was adjusted such that the Fe I and Fe II lines produced the same 
abundance within the errors. Microturbulence was determined by 
requiring that there was no significant slope in the Fe I abundance 
vs. predicted Fe I reduced equivalent widths based on the mean 
abundance of the lines. 

Table 2 lists SDSS ID, along with T eff , log(g) and [Fe/H] and 
their associated errors, plus microturbulent velocity and rotational 
velocity, where measurable, derived from our high-resolution spec- 
tra, and Table 3 lists SDSS ID and abundances of [a/Fe], [C/Fe], 
12 C/ 13 C, and loge(Li) determined from synthesis of the 6707A 
and 6104A lines. Alpha-element abundances are the error-weighted 
average of [Mg/Fe], [Ca/Fe] and [Ti/Fe], and the total error on 
[a/Fe] is the larger of the propagated errors or the dispersion in 
the three input abundances. 



4.2 Evolutionary states of the stars 

Previous searches for Li-rich stars have focused on g roups of stars 
with similarities: for example, Lebzel ter et al.l J2012I) selected only 
RGB stars in the Galactic bu lge, with metallicities in a narrow 
range, and lKumar et al.l feOllh use the Hipparcos catalog, in which 
the distances to all stars are known quite precisely, as their initial 
data set. However, our data set is heterogeneous in terms of age, 
composition and position in the Galaxy. Based on the metallicities 
and alpha-enhancements derived from the high-resolution spectra, 
we divide our Li-rich stars into four broad groups and assign rep- 
resentative ages. Figure 5 shows the four groups in the [a/Fe] 
- [Fe/H] plane, which is commonly used to compare the star- 
formation and chemical-enrichment histories of composite stel- 
lar popula tions in dwarf galaxies to the components of the Milky 
Way (e.g..lVenn et alj2003l : IShetrone et al.l2003l ; iKoch et alj2008l ; 
iKiribv et alj|201ll) . Briefly, moderate- and low-metallicity stars are 
assumed to be old, while higher-metallicity, non-a-enhanced stars 
are significantly younger. For the highest-metallicity stars in the 
sample ([Fe/H] > —0.2), we use isochrones to estimate the ages. 



Figure 6. Effective temperature and surface gravity from high-resolution 
analysis, with measurement errors, for our 23 Li-rich stars. Lower- 
metallicity stars are in the left panel, and higher-metallicity stars are in the 
right panel. The different symbols represent groups in evolutionary phase, 
and representative isochrones are overplotted. Horizontal lines mark the po- 
sition of the RGB bump on each isochrone, as calculated from our empir- 
ical relation given in Section 2.2. Triangles are RGB stars, circles are red 
clump stars, and squares are stars that have not yet reached first dredge-up 
(whose Li abundances are not relevant to the question of Li production on 
the RGB). The three plain crosses at the top of the left panel are stars whose 
log(g) and T e f f are not compatible with any reasonable isochrone, and we 
believe them to be post-AGB stars. 



Figure 6 demonstrates that our initial [a/Fe] - [Fe/H] age esti- 
mation was reasonable: both panels compare our values of T c ff and 
log(g) measured from high-resolution spectra against representa- 
tive isochrones, with older stars in the left panel and younger stars 
in the right panel, and the majority of stars lie near the giant branch 
corresponding to the appropriate initial values. The horizontal lines 
in each panel mark the log(g) at the RGB bump for each isochrone 
shown, calculated from our empirical calibration given in Section 
2.2. The squares represent stars that have not yet begun to ascend 
the RGB, and their high Li abundances are presumably primordial. 
The triangles represent RGB stars, both young and old, and the cir- 
cles show red clump stars, which all sit to the left of the isochrone 
at their metallicity, and are all at intermediate age. The three crosses 
at the top of the left panel are three low-metallicity stars with pa- 
rameters that are apparently inconsistent with the isochrones; we 
will address these stars in Sec. 5. The last two columns of Table 2 
give our final estimates for age and evolutionary phase. These val- 
ues are not given for the three stars with parameters incompatible 
with the isochrones. 

We are able to determine 12 C/ 13 C for all six of our red clump 
stars, using synthesis of the 4297A CH feature. In all of these stars, 
the ratio is fairly low ( 12 C/ 13 C < 15), indicating that first dredge- 
up has efficiently mixed the surface material with the stellar inte- 
rior. 

The Li-rich stars are mainly found low on the RGB, which 
might be expected from Li production associated with the onset 
of deep mixing. However, the presence of some bright Li-rich gi- 
ants and red clump stars indicates that it must also be possible for a 
star to produce Li during other evolutionary phases. Current models 
of thermohaline mixing on the red giant branch are fairly uncon- 
strained in terms of whether matter is transported away from the 
hydrogen-burning shell quickly enough for the Cameron-Fowler 
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Figure 7. Lithium abundance versus surface gravity for the stars in our 
sample. The typical value for loge(Li) decreases in the more evolved stars. 
Symbols are as in Figure 6, and there is no clear difference in the behavior 
of RGB stars and red clump stars in this plane. 



machanism to produce Li. Figures 7 and 8 show our resulting Li 
abundances versus log(g) and [Fe/H], respectively, using the same 
symbols as in Figure 6. 



5 DISCUSSION 

We have identified 23 new Li-rich post-main-sequence stars: 5 
bright red giants, 6 red clump stars, 6 RGB-bump stars, 3 lower- 
RGB stars and 3 post-AGB stars. Although schematic models exist 
for the production of Li at onset of deep mixing, near the RGB 
bump, there is not a clear theoretical explanation for the processes 
that produce and destroy Li in post-main-sequence stars. 

There are no targets in common between our data set and 
the six recent studies mentioned in the Introduction. However, 
our initial search of the DR7 CAS included SEGUE observations 
in Galactic globular a nd open clust e rs, all owin g a possible over- 
lap wi th the studies of Pila chowskil dl986h and Pilachowski et al. 
d 19881) . Our star SDSS 12356+5 633 is the known Li-rich giant NGC 
7789 K301 JPilachowskiH l986). It is chosen as a likely followup 
candidate according to our low-resolution selection criteria, and our 
high-resolution analysis confirms its high Li abundance. The fact 
that this star is in a cluster can be used to check our analysis. The 
mean metallicity of NGC 7789 is -0.04 ± 0.05 and 0.0 4 ± 0.07 
from iTautvaisiene et al.l d2005l) and IPancino et al.l ( |2010[) . respec- 
tively, which is c onsistent with our v alue of —0.05 ± 0.07. In addi- 
tion, for this star iPilachowskil dl986f) note that it is a weak G-band 
star, and we find that this star has one of the lowest carbon abun- 
dance ratios in our sample, [C/Fe]= —0.7. Further, the position of 
K301 on the CMD of NGC 7789 suggests that it may be an AGB 
star, or its low carbon abundance may cause the star to move bl ue- 
ward in the CMD via the Bond-Neff Effect dBond & Nefl 19691). 
We do not recover the Li-rich stars discussed in Roederer et al 



2008), Kumar & Reddv ( 2009), ICarnev et al.l lfl998). Kraft et al 



1999). lSmith et all Jl999h or lCharbonnel & Balachandranl fcOOO) 



which were not observed by SEGUE. Two of the Li-rich stars we 
selected for high-resolution followup are serendipitously located 
in the Kepler field. SDSS 11901+3808 is KIC 2968828, which 
has asteroseismology-based stellar parameters of log(g)= 2.393, 
T e ff = 4270, [Fe/H]= -0.823. While the surface gravity is a good 
match to our determination, the temperature and metallicity are 



Figure 8. Lithium abundance versus metallicity for the stars in our sample. 
The higher-metallicity stars show a broader range in loge(Li), and a higher 
mean, than the lower-metallicity stars. The rate of deep mixing is a sensitive 
function of stellar metallicity, so this may indicate that higher-metallicity 
stars are better able to produce lithium, or less able to destroy it. Symbols 
are as in Figure 6. 



quite different, and we believe that the discrepancy is likely due to 
uncertainties in the Kepler anaylsis techniques. SDSS J1909+3837 
is KIC 3531579, which does not have any publically available Ke- 
pler data yet. 

The three stars in our high-resolution data set that have stel- 
lar parameters incompatible with the isochrones are quite unusual. 
With their very low surface gravities, we identify them as post- 
AGB stars. The metallicities we derive for two of them are quite 
low ([Fe/H]~ —2.5). With these very low metallicities, these two 
stars are marginally consistent with being RGB-tip stars, but the ex- 
tremely low gravities are also consi stent with a known population 
of metal-deficient post-AGB stars dParthasarathvl fooO). We note 
that it has been hypothesized that very low metallicities in post- 
AGB stars are a result of dust formation in extremely cool atmo- 
spheres and are not representative of the initial stellar composition. 



5.1 Internal and external sources of lithium 

The initial photometric association of Li-rich stars with the 
RGB bump and the red clump led naturally to a suggestion 
dCharbonnel & Balach andran 2000) that the destruction of the /i- 
barrier created by first dredge-up causes a brief episode of Li pro- 
duction. This cannot be the explanation for high Li abundances in 
bright giants, far from the RGB bump, but any process that causes 
rapid mass motion between the envelope and the hydrogen-burning 
shell could in principle produce Li at the surface via the Cameron- 
Fowler mechanism: 3 He(a, 7) 7 Be(e~ ,i/) 7 Li. 

One potential mechanism for such stirring is proposed by 
iDenissenkovl d2012h . who note that, in their models of red giants 
with rapid interior rotation, the temperature gradient in the region 
between the hydrogen-burning shell and the envelope can become 
quite large, permitting deep mixing to operate much more quickly 
than usual. One consequence of this accelerated deep mixing is that 
the photometric loop executed during the RGB bump is much more 
extended, pushing RGB -bump stars into the same color-magnitude 
region as red clump stars. This indicates that perhaps the Li-rich 
stars that have been identified as red clump stars (including the six 
in this study) are actually RGB-bump stars with rapid internal ro- 
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Table 2. ID and stellar parameters from high-resolution analysis, for the 23 stars for which we could determine high-resolution stellar parameters 



SDSS ID 


T e ff 


17 


log(g) 


a 


[Fe/H] 


a 


v t urb(kms 1 ) 


a 


v rot (kms -1 ) 


Age (Gyr) 


EP 








O. u 






U. 1 1 








Q 


i ' 1 1 1 


SDSS 19^^-4-^798 

lJ-L/iJlJ JiJJJTJ / Z.O 


5025 


76 


3 


0.21 


+0.23 


0.08 


1.67 


0.11 


... 


4 


RC 


kJi_-/OiJ JUOWO uOU 


4900 


103 


3 


0.41 


+0.14 


10 


1.72 


0.14 




2 


RGB 


SDSS TO^m 4-71 


5075 


63 


3 1 


0.18 


+0.14 


06 


1.79 


0.08 




1 


RC 


iJi_-/iJiJ J 1 7TUU 1 Z 


4800 


100 


3 


0.41 


+0.06 


0.10 


1.70 


0.13 


... 


3 


RGB 


SDSS TOfW 4-4.0S9 


4900 


62 


2 9 


0.22 


+0.04 


0.07 


1.73 


0.08 


... 


2 


RC 


SDSS T074.S4-7109 

kJi_-/OiJ JLlZ.'+.JT/ 1\JZ, 


4850 


70 


2 1 


0.25 


+0.01 


07 


1.77 


0.07 




4 


RC 


SDSS T0f^74-76f)4. 


5050 


81 


2 8 


0.35 


+0.00 


09 


1.84 


0.11 




4 


RC 


SDSS T93Sfi-i- < !fi'3'3 

v3i_-/iJij JiJJUTJUJJ 


4750 


73 


2 5 


0.26 


—0.05 


0.07 


2.00 


0.07 


... 


j 


RGB 


iJi_-/iJiJ JUJJJTUJlt 


5200 


61 


2 8 


0.22 


—0.08 


0.06 


1.66 


0.08 




4 


RC 


SDSS J2200+4559 


4675 


100 


2.4 


0.37 


-0.12 


0.11 


1.99 


0.39 


6 


2 


RGB 


SDSS J0304+3823 


5125 


64 


2.6 


0.24 


-0.20 


0.07 


1.33 


0.09 




4 


RC 


SDSS J1901+3808 


4450 


86 


2.3 


0.45 


-0.29 


0.08 


1.78 


0.09 




4 


RGB 


SDSS J0720+3036 


5250 


142 


3.4 


0.52 


-0.29 


0.15 


1.40 


0.25 


10 


4 


SGB 


SDSS J1909+3837 


4450 


123 


2.2 


0.47 


-0.30 


0.13 


2.13 


0.13 




4 


RGB 


SDSS Jl 105+2850 


5625 


156 


3.9 


0.58 


-0.46 


0.18 


1.05 


0.31 




4 


SGB 


SDSS J0654+4200 


5200 


50 


2.6 


0.18 


-0.47 


0.05 


1.88 


0.06 


12.5 


4 


RC 


SDSS J1607+0447 


5175 


100 


2.0 


0.34 


-1.37 


0.13 


1.49 


0.15 




10 


RGB 


SDSS J1310-0012 


4550 


66 


1.0 


0.28 


-1.54 


0.06 


1.74 


0.07 




12 


RGB 


SDSS J0936+2935 


5125 


120 


0.0 


0.60 


-1.62 


0.15 


2.59 


0.23 








SDSS J1522+0655 


4975 


140 


1.9 


0.58 


-2.34 


0.15 


1.89 


0.18 




12 


RGB 


SDSS J0831+5402 


4725 


106 


0.2 


0.66 


-2.46 


0.12 


2.32 


0.23 








SDSS J1432+0814 


4625 


38 


0.3 


0.3 


-2.65 


0.06 


1.90 


0.13 


9.1 







Table 3. ID and abundances from high-resolution analysis, for the 23 stars for which we could determine high-resolution stellar parameters 



SDSS ID 


[a/Fe] 




a 


[C/Fe] 




a 


12 C/ 13 C 


loge(Li) 6708 


a 


loge(Li) 6104 




17 


SDSS J2206+4531 


-0 


29 





07 


-0.5 





3 




2.80 


0.30 


< 2.9 






SDSS J2353+5728 


-0 


08 





06 


-0.5 





2 


8 


3.10 


0.15 


< 3.2 






SDSS J0808-0815 


-0 


34 





18 


-0.6 





2 




2.45 


0.15 


< 2.7 






SDSS J0301+7159 


-0 


15 





08 


-0.6 





1 


15 


2.70 


0.10 


< 2.9 






SDSS J20 19+60 12 


-0 


08 





06 


-0.7 





2 




2.77 


0.20 


< 2.8 






SDSS J0652+4052 


-0 


05 





06 


-0.4 





2 


15 


3.30 


0.20 


3.0 





3 


SDSS J0245+7102 


-0 


07 





01 


-0.4 





1 


15 


2.58 


0.10 


2.4 





•1 


SDSS J0632+2604 


-0 


08 





05 


-0.4 





1 




4.20 


0.2 


4.1 





2 


SDSS J2356+5633 


+0 


06 





05 


-0.7 





1 




2.88 


0.15 


2.9 





2 


SDSS J0535+0514 


+0 


06 





04 


-0.5 





2 


15 


2.90 


0.10 


2.7 





3 


SDSS J2200+4559 


+0 


05 





08 


-0.8 





2 




3.05 


0.20 


3.0 





2 


SDSS J0304+3823 


-0 


01 





01 


-0.3 





2 




2.40 


0.10 


< 2.6 






SDSS J1901+3808 


+0 


17 





07 


-0.7 





1 




2.52 


0.15 


2.4 





3 


SDSS J0720+3036 


-0 


23 





13 


-0.4 





2 




4.55 


0.20 


4.2 





2 


SDSS J1909+3837 


-0 


20 





07 


< -0.4 








1.95 


0.15 


< 2.2 






SDSS Jl 105+2850 


+0 


09 





48 


-0.4 





4 




3.24 


0.30 


< 2.9 






SDSS J0654+4200 


+0 


02 





06 


-0.4 





1 


1 


3.30 


0.20 


3.1 





3 


SDSS J1607+0447 


+0 


33 





06 


-0.6 





2 




2.55 


0.15 


2.7 





3 


SDSS J1310-0012 


+0 


32 





09 


-0.7 





2 




2.15 


0.15 


2.2 





3 


SDSS J0936+2935 


+0 


38 





11 


0.3 





4 




2.70 


0.30 


< 3.0 






SDSS J1522+0655 


+0 


50 





20 


-0.7 





3 




2.25 


0.20 


2.6 





4 


SDSS J083 1+5402 


+0 


33 





11 


0.1 





3 




2.55 


0.20 


< 2.7 






SDSS J1432+0814 


+0 


12 





10 


-0.2 





2 




4.02 


0.30 


3.6 





3 



tation. Asteroseismic data on these stars could clarify their evolu- 
tionary phase. This process is an interesting possibility, but does 
not explain Li-rich bright RGB stars. 

The engulfemen t of a giant planet has also been suggested 
dCarlberg et ap2012l) as a way to "replenish" the surface Li abun- 
dance of a red giant. This process should be distinguishable from 
intrinsic Li production by deep mixing, since the 12 C/ 13 C ratio 
should rise following planet accretion, but should fall as a result of 



deep mixing. This is a process that could in principle happen at any 
point in stellar evolution, though the dramatic radial expansion that 
occurs as a star ascends the RGB would make it far more likely to 
happen during that phase. 
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6 SUMMARY 

The work presented here demonstrates the potential of large spec- 
troscopic surveys such as SEGUE as a source of unusual targets 
for followup observations. The ability to select our initial data set 
from a very broad parameter space allows us to search for Li-rich 
stars in a relatively un biased way. Other f ollowup studies based 
on SEGUE data (e.g., ICaffau et alj l201ll ; iBonifacio et al.l l201ll ; 
lAoki et alj|2012h follow the same approach: an initial search in a 
large, multidimensional parameter space, which introduces some 
heterogeneity in the data set, but can avoid biases inherent in more 

focused studies. 

Although we agree with lKirbvetail d2012h that Li-rich field 
giants are perfectly normal stars in most respects, we must disagree 
with t heir conclusion, originally expressed by Ide La Reza et al.l 
( 1996), that the onset of deep mixing produces a brief phase of 
Li enrichment in all stars. Such a universal process would imply 
that Li-rich stars should only be found near the RGB bump, which 
they are not. We propose instead that deep mixing is capable of pro- 
ducing Li enrichment at any point above the RGB bump, based on 
thermohaline mixing models. This does not explain what triggers a 
phase of Li production, but does fit comfortably with the the wide 
range of evolutionary phase that Li-rich stars occupy, and the short 
duration of the phase, which is inferred from the rarity of these 
stars. 

The combination of the recent large-sample studies of Li- 
rich stars is extremely important for advancing our understand- 
ing of the origins of these stars. Primarily, the new data demon- 
strates the ubiquity of Li-rich stars across Galactic and extragalac- 
tic environments - Li production is clearly a stochastic and short- 
lived, possibly recurrent event in stellar evolution. The next sig- 
ni ficant step will need t o take an approach similar to the study 
by ICarlberg et al.l d2012l) : identifying unique features of the var- 
ious proposed sources for Li enrichment in post-main-sequence 
stars, such as correlated variations in particular elemental or iso- 
topic abundances, and searching specifically for those signals. 
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APPENDIX A: COMPARISON OF HIGH-RESOLUTION 
AND SSPP STELLAR PARAMETERS 

During the course of this project, revised SSPP stellar parame- 
ters (including effective temperature, surface gravity and metallic- 
it y) were published for our stars as p art of SDSS Data Release 
9 dSDSS-III Collaboration et alj|201 2h . Since we have derived at- 
mospheric parameters from our high-resolution spectra, we present 
here a brief study of the comparison between the DR9 SSPP and 
high-resolution parameters, for the 20 stars with reasonable high- 
resolution parameters. High-resolution stellar parameters are listed 
in Table 2; DR9 SSPP parameters for the same stars are given in 
Table Al. 

Figures Al, A2 and A3 compare the high-resolution values for 
Teff, log(g) and [Fe/H], respectively, to the DR9 SSPP values, 
for the 20 stars with reasonable high-resolution parameters. There 
is scatter, but not a significant trend, between the high-resolution 
and SSPP determinations of T e ff and [Fe/H], but we find a sig- 
nificant trend with little scatter between high-resolution and SSPP 
determinations of log(g). We also find a significant trend, with a 
slope of 102±39 K/dex, between TgP -T^| pp and [Fe/H] SSPP , 
indicating that the SSPP algorithms are somewhat susceptible to 
temperature-metallicity degeneracy. The fact that we use different 
iron lines to determine [Fe/H] for the metal-poor and metal-rich 
stars in our sample could also contribute to this trend, but we would 
not expect that to cause an effect as large as what is observed. 
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Figure Al. Difference between high-resolution and DR9 SSPP determina- 
tions of T e ff, as a function of high-resolution T & ff. The best-fit dashed 
line is consistent with no offset between the two temperature scales, but 
significant scatter. Symbols are as in Figure 6. 
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Figure A2. Difference between high-resolution and DR9 SSPP determina- 
tions of log(g), as a function of high-resolution log(g). The best-fit dashed 
line showes a significant trend. The SSPP is known to have difficulty deter- 
mining log(g) in low-gravity stars, since the indicators it uses lose sensi- 
tivity; however, the trend to high gravity is unexpected. Symbols are as in 
Figure 6. 
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Table Al. ID and stellar parameters from the DR9 SSPP, for the 23 stars for which we could determine high-resolution stellar parameters 
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Figure A3. Difference between high-resolution and DR9 SSPP determina- 
tions of [Fe/H], as a function of high-resolution [Fe/H]. A best-fit line is 
shown; however, this distribution is consistent with a small offset between 
the two metallicity scales and ~ 0.3 dex scatter. Symbols are as in Figure 
6. 



